arXiv:nucl-th/0211032v2 4 Apr 2003 


Collective Properties of Low-lying Octupole 
Excitations in 20^ a 40 and l 8 0 20 

X.R. Zhou a,b , E.G. Zhao a,b,d , B.G. Dong c , X.Z. Zhang c , 

G.L. Long a,d 

a Department of Physics, Tsinghua University, Beijing 100084, P R- China 

b Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100080, 

P.R. China 

c China Institute of Atomic Energy, P.O. Box 275, Beijing 102431, P.R. China 

d Center of Theoretical Physics, National Laboratory of Heavy Ion Accelerator, 

Lanzhou 730000, P.R. China 


Abstract 

The octupole strengths of three nuclei: (3 —stable nucleus g 2 8 -P 6 i 26 ) neutron skin 
nucleus ^Ca^o and neutron drip line nucleus g 8 O 20 are studied by using the self- 
consistent Hartree-Fock calculation with the random phase approximation. The 
collective properties of low-lying excitations are analyzed by particle-vibration cou¬ 
pling. The results show that there is the coexistence of the collective excitations and 
the decoupled strong continuum strength near the threshold in the lowest isoscalar 
states in ^Ca^o and g S O 20 - For these three nuclei, both the low-lying isoscalar states 
and giant isoscalar resonance carry isovector strength. The ratio B(IV)/B(IS) is 
checked and it is found that, for g 2 8 P&i 26 , the ratio is equal to ( N f z ) 2 in good 
accuracy, while for 20^040 and g 8 O 20 > the ratios are much larger than ( N j A Z ) 2 - The 
study shows that the enhancement of the ratio is due to the excess neutrons that 
have small binding energies in 2oC a 40 and | 8 0 2 o- 
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1 Introduction 


Various exotic properties are expected for nuclei far from ( 3 — stability. The 
collective properties of neutron drip line nuclei are especially interesting, be¬ 
cause neutrons with small binding energies show a unique response to external 
fields. 
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In the Ref. [1-7], monopole, dipole and quadrupole isoscalar (IS) and isovector 
(IV) giant resonances in stable and drip line, particularly neutron drip line 
nuclei, were studied by using the self-consistent Hartree-Fock (HF)calculation 
plus random phase approximation (RPA) with Skyrme interaction, and both 
the IS and IV correlations were taken into account simultaneously. It was 
found that for both /^-stable and drip line nuclei the giant resonances can be 
well described by the collective model [8,9]. For neutron drip line nuclei, how¬ 
ever, there is appreciable amount of low-lying strengths just above threshold 
and these low-lying strengths are nearly pure neutron unperturbed excita¬ 
tions. For example, in Ref[4], the quadrupole strength of the neutron drip 
line nucleus g S 02 o was analyzed and it was pointed out that there exists a so- 
called threshold strength, which is not of collective character and comes from 
the excitations of excess neutrons with small binding energies. In a recent 
publication [11], the low-lying octupole excitation of the neutron skin nucleus 
2 oCV 40 was studied and it was pointed out the low-lying (A N = 1) IS octupole 
states appear as collective excitation and are shifted down to very low energy 
region, due to the disappearance of the N=50 magic number. Low-lying oc¬ 
tupole excitations usually consist of transitions from occupied states to bound 
states, resonance states and nonresonance states. In this paper, we make a 
detailed study of low-lying (AN = 1) octupole excitations of the (3 —stable 
nucleus g]] 8 Pb V2 % , the neutron skin nucleus and the drip line nucleus 

g 8 O 20 - If can be seen that in all these nuclei, the low-lying octupole states 
near threshold appear as the coexistence of the collective excitations and the 
decoupled strong continuum strength. The properties of these low-lying states 
can be understood from the point of view of particle-vibration coupling. 

This paper is organized as follows. The theoretical formalism of HF plus RPA 
calculation is described in section 2. Numerical results and discussions are 
shown in section 3. A summary and conclusions are given in section 4. 


2 Formalism 


The unperturbed strength function is defined by[12] 


So = J2\ W A |°> 1 2 S(E-Ei) 

i 

= -Im Tr(Q^G 0 (E)Q x ) 

7r 

while RPA strength function is given by 


(i) 


s = Y. I W<2» \ 2 S(E-E„) 

n 


2 



(2) 


= -ImTr(Q x ^G RPA (E)Q x ), 

^T 

where Go is the noninteracting p-h Green function, and Grpa(E) is the RPA 
response function including the effect of the coupling to the continuum, 

Grpa = Go + GoVphGRPA = (1 — Gov p h ) 1 Go- (3) 

In Eqs. (1) and (2), Q x represents one-body operators and is written as 


QT 3 ' r -° = En 3 I 3 for isoscalar octupole strength, (4) 

i 

QT 3 - r=1 = £ T z (i)rfY 3lM (ri ), for isovector octupole strength. (5) 

i 

The transition density for an exited state | n) is defined by 

Sp n0 (r) = {n\J2^r-fi)\0), (6) 

i 

which can be obtained by RPA response function 

Spno(r) = a j Im[GR PA (r, r* ; E res )}Q x (r > )dr* , (7) 


where the normalization factor a is determined from the transition strength 
£(A) by 


a = 


1 

Gm)' 


( 8 ) 


The radial transition density is defined by 
Sp„o(r) = Sp x (r)Y^(f). 


( 9 ) 


The transition current is 

Jno(r) = (n\ E)Vi ~ vAr- n)}|0), (10) 

which can be expanded in the complete set of the vector spherical harmonics, 

•MO = (—0 £ Jx,T M jr). (11) 

1 =Ail 
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The radial current component J\i is defined as 


J\i — i j t xl ^(f)J n0 (r)dr 

= < n IE - r *)[W*) x (v* + - v-)M I o). (12) 

In the Bohr-Tassie model, the transition density [8,9] is 

Sp\ T (r) oc r x ~ l — for A > 0, (13) 

dr 

and the radial current components [13] are 

! r x ~Vo( r ), for l = X — 

0, for l = A + 

The Tassie transition density (13) is 
relationship 

S( A) = | J 5p\ T (r)r ( ' X+2 ' 1 dr \ 2 , (15) 


(14) 

1 . 

normalized according to the following 


where S( A) is the transition strength of RPA state. This means that the nor¬ 
malized Tassie transition density should give the same strength as that of RPA 
state. 

In the present work, the properties of low-lying octupole excitations are stud¬ 
ied from particle-vibration coupling. For octupole excitation, we use the radial 
dependence of the particle-vibration coupling 

Fpv(r) ~ r i dU ( r \ (16) 

dr 

or 

Vp V (r) ~ (17) 

where U[r) is the HF potential and /?o( r ) is the ground state density. Eq.(16) 
has been successfully used for the coupling of particle to shape oscillations in 
Ref. [8], 
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The sign of the ratio 


< p\Vpv(r)\h > = J fip ph (r)Vpv(r)r 2 dr g 

< p|r 3 |/r > / 5p p h(r)r 5 dr 

determines the influence of particle-vibration coupling on the strength of the 
unperturbed p — h excitations[14], The magnitude of the ratio is a measure 
of how strongly the unperturbed strength of this p — h excitation is modified 
by performing RPA calculation. If the ratio is equal to zero, the unperturbed 
strength of this p—h transition will remain unchanged by the RPA correlation. 


3 Results and discussions 


We first perform the HF calculation with the SkM* interaction and then use 
the RPA with the same interaction including simultaneously both the IS and 
the IV correlation. It is solved in the coordinate space with the Green’s func¬ 
tion method, taking into account the continuum exactly. 

The calculated unperturbed strengths of g 2 8 P 6 i 26 , and g 8 02 o are shown 

in Fig. 1. In the harmonic oscillator model, the low-lying octupole strength 
(AN = 1) is approximately equal to the high-lying (AN = 3) strength in the 
large N limit, and the low-lying and high-lying octupole strengths approxi¬ 
mately exhaust 25% and 75% of the energy weighted sum rule, respectively. 
In real (3 —stable nuclei N is finite and the low-lying strength is less than that 
predicted by harmonic oscillator model in the large N limit. From Fig. 1(a) it 
can be seen that for |) 8 Pb 12 6 , the unperturbed low-lying (AN = 1) strength 
is centered at about 8 MeV and spread over the same order of energy range. 
It exhausts approximately 40% of total strength. The high-lying (AN = 3) 
strength is centered at about 24 MeV and exhausts approximately 60% of 
total strength. For 20 Ca 40 , the unperturbed strength (see Fig.l(b)) spreads 
over an energy range from 2 MeV to 15 MeV and exhausts about 60% of 
total strength. While for | 8 0 2 0 (see Fig.1(c)), except a few high-lying proton 
excitations, nearly all of the octupole strength lies within the low energy re¬ 
gion. For 20 G(Z 40 and large amount of neutron unperturbed octupole 

strengths are shifted down to very low energy region. This downward shifting 
of the octupole strengths is attributed to the disappearance of magic number 
N=50 for 2 q CCI 4 .Q as shown in Fig.2(a), and the disappearance of magic num¬ 
bers N=20 and N=28 for | 8 02 o as shown in Fig.2(b). For all the three nuclei, 
the calculated energy weighted sum rules are equal to the classical sum rules 
to a good accuracy. 

In Fig.2, one-particle energies are given with fixed neutron number for (a)N=20 
and (b)N=40 as a function of proton number. It is easy to see that in Fig.2(a), 
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when it is near the neutron drip line, the magic number N=20 is disappearing, 
which agrees with the experimental observation[15], and new magic number 
N=16 appears, which is consistent with the conclusion in Ref. [16]. Similar 
phenomena appear also in Fig.2(b). Near the neutron drip line, N=50 magic 
number is disappearing. 

Fig .3 gives the IS and IV RPA octupolc strength for the gSfP&m, 2 o^ a ^o 
and g 8 0 2 o- For (3— stable nucleus g 2 8 F > ^i 26 (Fig.3(a)), we see a strong IS peak 
below threshold at Ex=3.48 MeV with strength B{ A = 3, IS : 3 - —> 0 + ) = 
5.68 x 10 5 /m 6 , which corresponds to B(E3 : 3~ —> 0 + ) = (^f ) 2 x 5.68 x 
10 5 /m 6 = 0.89 x 10 5 e 2 /m 6 . These values are comparable to experimental 
data B(E3 : 3“ —■> 0 + ) = 1.0 x 10 5 e 2 /m 6 at energy 2.61 MeV . The IS 
giant resonance is at about 20.5 MeV and the IV giant resonances are mainly 
distributed in the energy region from 25 MeV to 35 MeV. 

In Fig.4 and Fig.5 we show,for the case of the transition densities and 

radial current components of IS excitation at Ex=3.48 MeV, IS giant resonance 
at Ex=20.5 MeV, and IV giant resonance at Ex=34.2 MeV, together with the 
prediction of the Bohr-Tassie model, respectively. From these figures, it can be 
seen that the IS and IV octupole excitations are surface modes and they are 
well described by the collective model. For the strong collective IS state at 3.48 
MeV in Fig.4(a) and Fig.5(a), there are some difference between two models 
at small but in the surface region with r at about 7 fm, this excitation can 
well be described by the Bohr-Tassie model. 

In neutron-excess nuclei, the IS mode (shape vibration) gives rise to a IV 
moment which is proportional to ( N — Z), so the strength of an IS mode 
carries IV strength, and the ratio of IV strength/IS strength is expected to 
be (^--p -) 2 as pointed in Refs. [4,10,17]. We calculated this ratio for collective 
IS modes at Ex=3.48 MeV and giant IS resonance around Ex=20.5 MeV in 
g 2 8 P 6 i 26 - Both modes give the ratios (^§yp ) 2 to a good accuracy. 

For 20 CV 40 (Fig.3(b)), the main IS strengths are shifted down to low energy 
region. Below 5 MeV there are 4 strong IS peaks and one of them is below 
threshold (Ex=1.91 MeV). ^Frorn the transition densities and radial current 
components in Fig. 6 , the IS collective state at Ex=1.91 MeV can be described 
by the Bohr-Tassie model. The transition densities and radial current compo¬ 
nents of the other three IS peaks below 5 MeV are shown in Fig.7 and Fig. 8 , 
respectively. 4 From Fig.7 it can be seen that these three IS peaks are surface 
modes and both neutrons and protons contribute. This indicates they are col¬ 
lective. In the central region of the nucleus there are differences between the 
transition densities of the Bohr-Tassie model and those of RPA calculation. 
But in the surface region the results of RPA calculation are similar to those 
of collective model. From Fig .8 we see that for the currents, the Bohr-Tassie 
model results differ from the ones from RPA substantially, even in the surface 
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region. For example, in our RPA calculation the small component . 73,4 (r) is 
comparable to the large component j 3 j 2 (r), but the Bohr-Tassi model gives 
_ 7 3 ) 4 (r) — 0. This difference will be explained from the viewpoint of particle- 
vibration coupling later. 

The ratios of IV strength/IS strength for these 4 strong IS peaks (below 5 
MeV in 20 Ca 40 are also calculated. We find that, however, unlike the case 
in (3— stable nucleus Pb i 26 , the calculated ratios are larger than (^A ) 2 
by factor 2 to 4. This result is closely related to the neutron orbits of small 
binding energies and they push down the unperturbed octupolc strengths to 
very low energy region. 

For | 8 0 2 0 (see Fig.3(c)), the high-lying IS giant resonance nearly disappears 
and almost all the IS octupole strengths are shifted down to low energy. The 
transition densities and radial current components for the lowest three peaks 
are shown in Fig.9 and Fig. 10, respectively. From Fig.9 we see that the three 
peaks are surface modes and both protons and neutrons contribute. In Fig.l 
and Fig.3 we notice that the IS strengths of these peaks are larger than those 
of unperturbed ones, which shows they are collective. The Bohr-Tassie model 
can approximately describe the transition densities in the surface region. The 
Bohr-Tassie model gives quite different currents properties for these IS peaks 
in Fig.10. The reason will be analyzed later. 

We also calculated the ratios of IV strength/IS strength for these three IS 
peaks. The calculated values are factor 3 to 5 larger than (^2 ±) 2 j n | 8 O 20 . 
Similar to ^Ca^o, this is also related to the least bound neutrons. 

Next we try to understand the collective properties of these low-lying IS exci¬ 
tations in 2 oCa 4 o and § 8 0 2 o based on the particle-vibration coupling. We have 
checked that the signs of the ratio Eq. (18) in our cases are always positive, so 
we only consider their absolute values here. Fig.11 shows a few low-lying neu¬ 
tron unperturbed octupole strengths in 2 QC'a 4 o for radial operators r 3 , r 2 dU ^ 
and r 2 dp0 ^ , and Fig.12 shows the corresponding quantities for g 8 0 2 o- Here 
where U(r) is the neutron radial HF potential and po( r ) is the ground state 
density. 7 From Fig. 11 we see that for the transition from bound state to non¬ 
resonance state, I /5 —» 3s 1 , there is a pronounced peak in octupole strength 
2 2 

function for radial operator ?’ 3 , but almost no peaks for the other two radial op¬ 
erators r 2 dU ^ and r 2 dpo d ^ , respectively. Just like the transitions from bound 
states to bound states, If 5 —> lgs and 2 p 3 —> lc/a, there are strong peaks of 
octupolc strengths for all three radial operators in corresponding energy re¬ 
gion for the transitions from bound states to resonance states, I /5 —> 2d 5 and 

2 2 

2pi —» 2d 5 , . It means that the ratios Eq. (18) for transitions to nonresonance 
states are much smaller than those for the transitions to resonance states and 
to bound states. In RPA calculation the unperturbed octupole strength of the 
radial operator r 3 for the transition from bound state to nonresonance state 
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If 5 —> 3si will hardly be affected, but the strengths of the radial operator 
2 2 

r 3 for the transitions from bound states to resonance states and to bound 
states will be strongly absorbed into the collective excitations. Because the 
unperturbed octupolc strength of the transition from bound state to nonres¬ 
onance state I/s —> 3si for radial operator r 3 is mainly distributed within 
.22 

the energy region from 3.0 MeV to 4.5 MeV, there is the coexistence of the 
collective excitations and the decoupled strong continuum strength near the 
threshold in these three lowest IS peaks. That is why the Bohr-Tassie model 
can only well describe the transition density of RPA calculation in the surface 
region for the low-lying excitations in but it can not well describe the 

calculated transition current. 

^From Fig. 12 we see similar results for g S 02o- For transition from bound state 
to nonresonance state, Id 3 —> 2p3, the ratio in Eq. (18) is much smaller than 
those for transitions from tound states to resonance states, 1(4/2 —> I/7/2 and 
2s 1/2 —>► 1 / 7 / 2 - In these three lowest IS peaks, there is the coexistence of the 
collective excitations and the decoupled strong continuum strength near the 
threshold. That is the reason why the Bohr-Tassie model can only approxi¬ 
mately describe the calculated transition density in the surface region for the 
three lowest IS peaks in g 8 0 2 o> but it can not well describe the transition 
current of our calculation. 


4 Summary and Conclusions 


The octupolc vibrations for the (3 —stable nucleus g 2 8 P&i 26 , the neutron skin 
nucleus and the drip line nucleus g 8 0 2 o are studied. It is found that 

the lowest IS excitation below threshold for the nuclei 82 8 -P6i26 and 2oCa4o, 
and IS and IV giant resonances of the /3-stable nucleus g 2 8 -P&i 26 can be well 
described by collective model, at least in the surface region. 

For the neutron skin nucleus 2 oCa 40 and the neutron drip line nucleus | 8 0 2 o, 
there exist strengths of transitions from bound states to bound states, reso¬ 
nance states, and nonresonance states in the low-lying unperturbed neutron 
octupolc strength (AN = 1). The strengths of transitions to nonresonance 
states are nearly unaffected and the strengths of other transitions are strongly 
absorbed into collective states by taking into account the RPA correlation. So 
there is the coexistence of the collective excitations and the decoupled strong 
continuum strength near the threshold in the lowest IS states. 

We also find that, for the (3 —stable nucleus glfPfe^e- both low-lying IS states 
and giant IS resonances carry IV component and the ratios of IV strength/IS 
strength are equal to to good accuracy, but for the neutron skin nu¬ 

cleus 20^040 and the neutron drip line nucleus g 8 02 o, these ratios for a few 



lowest strong IS excitations are much larger than ( N ~\ Z ) 2 ■ These results are 
closely related to the small binding energies of neutron orbits in these nuclei. 
The octupolc transitions from these orbits are mainly distributed in low en¬ 
ergy region, so the contribution to low-lying (A N = 1) octupole states from 
neutrons are much larger than those from protons. 

This work is supported by the National Natural Science Foundation of China 
under contact 10047001 and the Major State Basic Research Development 
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and H. Sagawa for providing us with continuum RPA program. 


References 

[1] I. Hamamoto, H. Sagawa and X. Z. Zhang, Phys. Rev. C 53 (1996) 765. 

[2] I. Hamamoto and H. Sagawa, Phys. Rev. C 53 (1996) R1492. 

[3] I. Hamamoto and H. Sagawa, Phys. Rev. C 54 (1996) 2369. 

[4] I. Hamamoto, H. Sagawa and X. Z. Zhang, Phys. Rev. C 55 (1997) 2361; J. 
Phys. G 24 (1998) 1417. 

[5] I. Hamamoto, H. Sagawa and X. Z. Zhang, Nucl. Phys. A 626 (1997) 669. 

[6] I. Hamamoto, H. Sagawa and X. Z. Zhang, Phys. Rev. C 56 (1997) 3121. 

[7] I. Hamamoto, H. Sagawa and X. Z. Zhang, Phys. Rev. C 57 (1998) R1064. 

[8] A. Bohr and B. R. Mottelson, Nuclear Structure, Vol. II (Benjamin, New York, 
1975). 

[9] L. T. Tassie, Australian, J. Phys. 9 (1956) 407. 

[10] F. Catara, E. G. Lanza, M. A. Nagarajan, and A. Vitturi, Nucl. Phys. A 614 
(1997) 86. 

[11] I. Hamamoto, H. Sagawa and X. Z. Zhang, Phys. Rev. C 64 (2001) 024313. 

[12] I. Hamamoto, H. Sagawa and X. Z. Zhang, Nucl. Phys. A 648 (1999) 203. 

[13] T. Suzuki and D. J. Rowe, Nucl. Phys. A 286 (1977) 307. 

[14] I. Hamamoto nad X. Z. Zhang, Phys. Rev. C 58 (1998) 3388. 

[15] D. Guillemaud-Mueller, et al, Nucl. Phys. A 426 (1984) 37; T.Motobayashi, et 
ah, Phys. Lett. B 346 (1995) 9. 

[16] A.Ozawa, T. Kobayashi, Y. Suzuki, K. Yashida and I. Tanihata, Phys. Rev. 
Lett. 84 (2000) 5439. 

[17] H. Sagawa, Phys. Rev. C 64 (2002) 064314. 


9 




Fig. 1. Unperturbed octupole strengths of (a) 


l° 2 8 Pb 12 6, (b) fjCa 40 and (c)| 8 O 20 . 
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Fig. 2. One-particle energies vary with proton number for (a) neutron number 
N = 20 and (b) neutron number N = 40. 
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Fig. 3. Isoscalar and isovector octupole strengths of (a) g 2 8 P 6 12 6 , (b) 2 oC a 40 and 
(c) | 8 0 2 o- 




Fig. 4. (left) Radial transition densities of low-lying and high-lying IS and IV oc¬ 
tupole modes of § 2 8 P&i 26 - 

Fig. 5. (right) Radial transition current components of low-lying and high-lying IS 
and IV octupole modes of g 2 8 R^i 26 - 


11 
























Fig. 6. Radial transition densities and current components of low-lying IS octupole 
modes of 2oC a 40- 




Fig. 7. (left) Radial transition densities of low-lying IS octupole modes of 2oC a 40- 
Fig. 8. (right) Radial current components of low-lying IS octupole modes of 2oC a 40- 
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Fig. 9. (left) Radial transition densities of low-lying IS octupole modes of g 8 O 20 - 
Fig. 10. (right) Radial current components of low-lying IS octupole modes of | 8 O 20 - 
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Fig. 11. Few low-lying unperturbed octupole strength in for radial operators 

(a) r 3 , (b) r 2 dU ^ and (c) r 2 dp< ^^ , where U(r) is the neutron radial HF potential 
and po(r) is the ground state density. 
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Fig. 12. Few low-lying unperturbed octupole strength in g 8 O 20 for radial operators 
(a) r 3 , (b) r 2 dU ( ^ and (c)r 2 ^^, where £/(?’) is the neutron radial HF potential 
and po(r) is the ground state density. 
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